CO 



CN 



o 



Tunneling Gap as Evidence for Time-Reversal Symmetry 
Breaking at Surfaces of High-Temperature Superconductors 



R. B. Laughlii 

Department of Physics, Stanford University, Stanford, CA 94305 

and 
q^ ■ Lawrence Livermore National Laboratory, P. 0. Box 808, Livermore, CA 94-550 



Abstract 



J> ■ It is argued that recent Josephson junction and point-contact tunneling exper- 



iments, interpreted as intended by their authors, indicate that time-reversal 



symmetry breaking occurs at surfaces of cuprate superconductors. The varia- 

Tt; ' ' ' 

ON ■ tion among experiments and the failure of previous searches to find T- violation 



are ascribed to disorder and effects of 3-dimensionality. The "anyon" approach 



6 

i-rt ■ to the t-J model is shown to predict a conventional BCS order parameter of 

d x 2_ y 2 +ie d xy symmetry, with e roughly 3 times the doping fraction 5, which 



is consistent with these experiments but not demonstrated by them. 
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The purpose of this Letter is to point out that time- reversal symmetry breaking jl], the 
key prediction of the "anyon" approach to the high-temperature superconductivity problem, 
may already have been demonstrated in a series of recent experiments conducted for other 
purposes. While it is wise to be cautious, particularly in light of previous failures to detect 
T-violation and the possibility that one or more of these experiments may later prove to 
be wrong or misinterpreted, the implication of the experiments as they now stand is clear. 
The relevant experiments are the photoemission || and light scattering j|] measurements 
of the gap anisotropy, the microwave measurements of the low-temperature conductivity 
||, the Josephson phase coherence experiments of Wollman et al. ||, and the scanning 
tunneling microscope measurements of Hasegawa et al. [0] reporting an intrinsic energy gap. 
The essence of the argument is that all of these except the last point to the occurrence 
of a conventional BCS order parameter of d x 2_ y 2 symmetry, the behavior commonly found 
in theories based on incipient antiferromagnetism of the conducting electrons ||. This, 
however, is fundamentally incompatible with the last experiment unless the order parameter 
is complex, which is impossible unless the ground state violates T. Thus the simultaneous 
occurrence of <i-wave superconductivity and an energy gap in any part of the sample, if true, 
constitutes definitive evidence for T-violation. 

A d x 2_ y 2 order parameter A& is distinguished from a conventional s-wave order parameter 
by sign reversals. In high-Tc superconductors this is most easily discussed in terms of the 
idealized electron energy band 

E° k = -2t [cos(Jfc*&) + cos{k y b)] (1) 

where to is an energy parameter and b is the bond length of a 2-dimensional square lattice. 
One imagines creating a fermi sea by filling the states with E^ < 0. Allowing the electrons 
to interact through weak near-neighbor spin exchange || then leads to a superconducting 
state for which the quasiparticles have energies 



E k = ±y/{E$)* + \A k \* (2) 

where 



A fe = A [cos^fe) - cos(k y b)] . (3) 

The gap parameter A& is positive in the x-direction, negative in the ^-direction, and zero 
at the nodes in between. In an s-wave superconductor A*, is the constant A . 

Of the experiments listed above, only the Josephson experiment of Wollman et al. can 
directly sense the sign reversal of A&, even in principle. It accomplishes this by detecting 
electric currents spontaneously generated in a loop of Pb wire connected between the x— and 
y— faces of a high-Tc superconductor. Although the result reported by Wollman et al. H 
is positive and consistent with the large body of circumstantial evidence for (i-wave pairing, 
it has not yet been reproduced and is quite controversial. In what follows we shall assume 
that this experiment is right. It must be emphasized that it is the only direct evidence we 
have for sign reversal of A&. 

This experiment is supported by considerable indirect evidence for <i-wave superconduc- 
tivity with nodes at A; = (±7r/26, ±7r/2fe), in particular by the presence in all samples of 
low-energy excitations in the superconducting state. The evidence is, unfortunately, com- 
plicated by the disorder effects that plague these materials. For example, Giaver tunneling 
finds states in the gap so commonly that the claim of Hasegawa et. al. [|7j to have observed 
a clean gap, the experiment motivating this Letter, is widely questioned. At the same time, 
the large zero-bias conductance seen in most tunneling experiments is commonly dismissed 
as an artifact of disorder at the tunnel contact ||. The T 2 deviation of the penetration 
depth from its zero-temperature value, and its crossover to linear— T behavior above about 
5 K, seen in high-quality samples [H| is consistent with d x 2_ y 2 superconductivity only if 
disorder that is difficult or impossible to eliminate from the samples |ll|] is assumed to exist. 



The same is true of the zero-temperature microwave conductivity |5|JT0||. The intrinsic heat 
capacity below 10 K is not known for any high-T c superconductor because the signal is al- 
ways swamped by a large Schottkey-like heat capacity |12| similar to that expected of a spin 
glass. That at least some of these excitations are intrinsic and attributable to a ci-wave node 
is indicated by several less accurate or model-dependent experiments, the most accessible of 



which is the angle-resolved photoemission work of Shen et al. || . This reports measureable 
changes to the quasiparticle energies Ef. resulting from cooling the sample through its su- 
perconducting transition except near this special value of k. The complex temperature and 
polarization dependence of inelastic light scattering is accounted for quantitatively by d x 2_ y 2 
superconductivity if reasonable assumptions are made about the relevant matrix elements 4 . 
The same is true for the extensive magnetic resonance data ||13|| . While the situation is 



still confusing and controversial, it is clear that a large body of experimental evidence is 
consistent with the simultaneous presence in all samples of both disorder and a <i-wave 
node. 

Let us now consider the tunneling experiment of Hasegawa et al. 0] . It is very important 
for our argument that this experiment was performed with a scanning tunneling microscope 
tip, that it reported a well-developed energy gap only at certain places on the sample, 
and that it was inconsistent with conventional thin-film tunnel junction experiments, which 
never reveal a full gap ||. Barring the possibility that it is an experimental artifact, such 
as coulomb blockade [Il| or an effect of anisotropic tunneling, this result implies that a 



genuine energy gap in the quasiparticle spectrum develops at some places on the surface. 
The cleanliness of the reported gap suggests that this gap is intrinsic. Why it should develop 
only on islands is an open question. One possiblity is that it is destroyed almost everywhere 
on the surface by disorder. That disorder has the potential to do this, provided that the 
order parameter has sign reversals, is well known |11|| . This explanation probably cannot 
account for the absence of the gap in the bulk, although further experimental studies are 
required to tell for certain. A more likely explanation in this case is that the gap is suppressed 
because of the magnetic fields it would create. However, regardless of the mechanism by 
which the gap is destroyed, its existence anywhere in the sample indicates T- violation, since 
an energy gap is forbidden as a matter of principle in a <i-wave superconductor with a real 
order parameter. 

Let us now consider the specific complex order parameter 



A chirai = A o [ cos (^ & ) _ 008(^6) + ie sin{k x b) sm{k y b)} (4) 

where e is a number, implicit in the anyon technique. Because the imaginary part of this 
order parameter has sign reversals, it is susceptible to destruction by non-magnetic disorder 



11 and thus consistent with the idea that e might be zero everywhere except on islands. 



That there should be such an order parameter was first suggested by Rokhsar |15 |, who 
also pointed out that it was potentially inconsistent with the known phase diagram of these 
materials. Since the real and imaginary parts of Al hiral lie in different irreducible represen- 
tations of the lattice point group, they cannot mix in the Ginzberg-Landau functional, and 
thus must develop as separate order parameters as the temperature is lowered. This would 
give 2 phase transitions, rather than the observed 1. However, the suppression of the the 
smaller d xy order parameter almost everywhere in the sample would resolve this paradox 
by preventing the acquisition of long-range order and killing the lower phase transition. 
It would also kill all effects of global T-violation, and account for the failure of previous 
searches to find them ||. 

Let us now show that this order parameter is implicit in the anyon approach to the t — J 
Hamiltonian 

J 



Ti-t-j = 

<j,k> 



E |-*£sW+9<V^} (5) 



and estimate its magnitude. It should be noted such an order parameter is not uniquely 
obtained this way, but also arises in certain simple BCS Hamiltonians |l(J. As usual, Sj 



denotes the spin operator for the j th site, < j,k > denotes a sum on near-neighbor pairs, 
with each pair counted twice, and J = 0.1 eV and t = 0.5 eV are the spin-exchange and 
electron hopping matrix elements. We assume that the lattice has N sites and M = N5 
holes. 

The anyon superconducting state for the t — J model is an electron wavefunction of the 
form 

|*>= E a lu ... /M \^...J M ) (6) 



where | £i,...,£m) is a basis wave function describing "holons" at sites £i,...,£m- This is 
given explicitly by |T7| 

N M 



-1, ■■•■, Z M 



) = TJ(1 - n nnji ) [](1 - n^ - n ial )\ ^ flux ) (7) 



where nj ff = CjCr c ja is the number operator for an electron of spin a at site j, and | ^fi ux ) 
is the ground state of the commensurate flux Hamiltonian 

Wflux = -t Y. Y. ex P \ i 9 i 1 - 8 )( x 3 ~ x k)(Vj + Vk)/b 2 \c) a c ku (8) 

<j,k> a ^ A > 

with N — M electrons. Tifi ux violates T, as is required for the holon basis to be defined 
17] . It is also possible to construct the basis using 7i*f lux , in which case T is violated in 
the opposite sense. The holons defined by Eq. (7) obey 1/2 fractional statistics [|l7]] in that 
varying the expansion coefficients in Eq. (6) to minimize the expected energy (\l/ \H t ~j\ \&) 
is equivalent, when 5 is small, to solving the fermion problem [[TJ|] 

h 2 e 

j 

Ihc^ zx (r 3 - r k ) 
A J~2~2^ \r--rA 2 [ } 

z e kjij \' 3 '*l 

with isospin. The latter corresponds to the valley degeneracy of the holon band structure, 
which in the 5^0 limit is described by Eq. (8) without spin [[T7J . Since the holon dispersion 
relation in this limit is given by 



^holon ^ ±2t^ C OS 2 (k x b) + COS 2 (k y b) , (11) 

we find that the valley minima occur at k — (0, 0) and (ir/b, 0) and are characterized by the 
mass m* ^ h 2 /{\/2b 2 t). 

The first step in the order parameter calculation is to compute the non-local order pa- 
rameter of the continuum anyon gas described by Eq. (9) ||18| . We will take the ground 
state to be 

I Canyon) = II { 3 ^^ WV*PP q P-q) } I ®Hf) (12) 



where p = 5/b 2 is the particle density, p q = £,,exp(ig • Tj) is the density operator, 
and | $>hf) is the ground state of the Hamiltonian obtained by substituting the mean-field 
vector potential < A > = hc/(2e)p yx for Aj in Eq. (9). The prefactor in this expression 
is the usual modification of the Hartree-Fock ground state implicit in the random phase 
approximation. Let if>+(z) and ip-{z) denote the operators annihilating holons with "up" 
and "down" isospin, respectively, at z = x + iy in the fermi representation. The nonlocal 
order parameter is given in terms of these by [|18| 

(Canyon IV'+^OV'-^OiV'+WV'+W + ^-{z)^-{z) \lp + (z 2 )lp-{z 2 ) | $anyon) 

(zl - z *) ( Z2 - z ) s . , 2 



-d 2 z 



[i.9 P y 



(13) 



| Zi — Z | | Z 2 — Z I I Z\ — z 2 1 — *■ oo 

The numerical value of 1.9 p reported here for the first time is obtained using the hypernetted 
chain technique fll7| . We will adopt the notation (if; + (z)if)-(z)) = 1.9 p as shorthand for this 
result. 

The second step is to convert this continuum order parameter to the site basis. The 
unitary transformation relating ip+{z) and ip-(z) to if} (J), the fermi operator to annihilate 
a holon at site rj = (£j,rrij)b, is simply the the matrix of 1-body eigenstates of Eq. (8) at 
the two valley minima. We thus have 



Hj) = v 7 ^ 



cos 



4 1 2N ' 

m„-)7r 

v 8 2 jJ 



ip+(z) + (-1)4 sin 



A 1 2, ' 

y 8 2 J ' . 



iP-(z) 



(14) 



In obtaining this expression, we have imagined the sample to be divided into 4-site cells, 
and that z defines the center of the cell containing j. Then, substituting Eq. (13) into Eq. 
(12), and using the fact that ip+{z) and if>-(z) anticommute with themselves and each other, 
we obtain 



(0 | tf (l)tf(2) | 0) ^ 1.9 ( J(_l)<i+ma+«i+<a)mi 

1 



v/2 V 



:-i) 



£i+£2+mi+m2 



■y _ f_n('i+«2)(mi+m2) A 



(15) 



The third step is to evaluate the matrix element of cLcj^ to annihilate a holon pair into 



the vacuum. This matter has been studied extensively in previous papers [17 and is too 
involved to discuss in detail here. We shall simply quote the result 



t J 



(£[,.. .,£' M Icfc^CjVfcJ^J £i,...,£ M ) 



M 



^e 1 j^e 2 k^e' 1 j'^e' 2 k / 11 ^ 



a=3 



(z'* + z'* k )/2-z* a ( Zj+Zk )/ 2 -z a 
|(4 + 4)/2 - * a | |(^ + * fc )/2 - z a \ 



,t J 



x (j 1 , k' \cf^Ck'i\ 0)(0 |c]|C{,J j, fc) + cyclic permutations 



(16) 



where | 0) denotes the state with no holons and Zj = (£j + irrij)b. The product on a in 
this expression is the same "unwinding" factor appearing in Eq. (12) and is the microscopic 
justification for its inclusion in Eq. (12). The remaining factor is given approximately by 



(o i4 T 4j 3 ,k) - \ Zl ^) .( <1 - €2)(mi+m2 ) ( _ 1) € 



i+mi 



X 



Z\ ~ Z 2 \ 

l-5)exp{-7r/4(l-<5)} 

<5exp {— -irS} 




1 and 2 near neighbors 
1 and 2 second neighbors 
otherwise 



(17) 



The final step is to combine Eq. (15) and (16) with Eq. (14) to obtain an expression 
for (c]-|-c^). This is accomplished by multiplying together the right sides of Eqs. (14) and 
(16) and then dividing out the factor (z\ — Z2)/\zi — z 2 \. This latter step accounts for 
the transformation of the basis functions defined by Eq. (7), which are symmetric under 
interchange of the holon positions, to their fermi representation. We obtain finally 



(c\^ 2i ) ^ j(*-^)("*-™i)(-l) 



m\— rr%2 



x 1.95 



\/2(l — 5) exp {— 7r/4(l — 5)} ; 1 and 2 near neighbors 

— 25 exp {— ti5} ; 1 and 2 second neighbors 

; otherwise 



This is equivalent to Eq. (4) with 



V85 (1 



(1-55)tt/4 



1-5 



(19) 



Let us now make some comments about this result. The first is that it is quite crude 
and should be compared with experiment carefully. For example, it makes the unphysical 
prediction that superconductivity occurs in the t-J model at any value of J, t, and 5. A more 



8 



formal development of the anyon approach using gauge theory techniques |19j cures problems 



of this kind, but is too technical to discuss here. Let us simply state the main results, i.) For 
realistic values of J/t the order parameter is significantly smaller than Eq. (17) predicts, due 
to retardation effects, ii.) The equations have an antiferromagentic instability for 5 < .05 
that depends weakly on J/t. iii.) The commensurate flux band structure of Eq. (6) has a 
gap collapse at S = 1/3 which causes calculations for doping fractions larger than this to 
be unreliable. That the calculation finds conventional superconducting order and <i-wave 
symmetry is not surprising. The anyon approach is a legitimate variational technique for 
the t-J model, which is known by more reliable methods to have a tendency to d-wave 



pairing ||20| . Also, it has been known for several years that "flux" vacua are fundamentally 
related to <i-wave superconducting states plf . The calculation is significant mainly because 
it predicts that superconducting pairing by means of "spin fluctuations" tends naturally to 
an order parameter with d x 2_ y 2 + ie d xy symmetry, with e significantly large. 
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